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We discuss the effect of dynamical charged particles in ^/(l)^ -1 Abelian effective theories of 
QCD. Screening and confinement at zero temperature are explained qualitatively. The effect of 
dynamically charged particle is expressed in terms of an effective action of integer electric currents. 
When a screening is expected, the static potential shows a flattening in the long-range region and 
a linear behavior in the intermediate region. We show why the screening is better observed in the 
Polyakov-loop correlators than in the Wilson loops. The breaking of the adjoint string is explained 
without the Z(N) picture. 
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Understanding of non-perturbative aspects of strongly 
coupled field theories - such as QCD - is always a chal- 
lenging task. From the study of confinement in QCD, it 
has been known that duality and topology are two impor- 
tant keywords. Namely 't Hooft idea Q of monopole con- 
densation after Abelian projection to U (l)^ 1 of SU (N) 
QCD has been shown to be successful in the infrared 
region. Abelian and monopole components seem to be 
dominant in the long-range region of quenched QCD §. 
An infrared effective monopole action has been derived in 
the continuum limit after a block-spin transformation of 
monopole currents || Q]. It is a quantum perfect action 
described by monopole currents. The monopole action 
shows that monopole entropy dominates over the energy. 
That is, monopole condensation occurs in the confine- 
ment phase. Using a quantum perfect operator we 
can evaluate the string tension analytically and prove 
the continuum rotational invariance of the static poten- 
tial |. 

Although the 't Hooft scenario seems correct in many 
respects, there remains an unsolved serious problem 
even in quenched QCD and also in full QCD. It is the 
screening-confinement problem extensively discussed in 
recent publications JtJ. How can screening be explained 
in this scheme? Consider a static pair of charge 2 par- 
ticles in quenched QCD. Since off-diagonal gluons can 
screen the charge 2 particle, the static potential is ex- 
pected to show flattening in the long-range region, i.e., 
screening. Actually, an adjoint Wilson loop and a charge 
2 Wilson loop show such behaviors in the strong coupling 
region. This problem was discussed extensively in Ref . |8) . 
The theory in terms of Abelian link fields or Abelian 
monopole currents alone becomes highly non-local if we 
integrate out off-diagonal gluon fields after an Abelian 
projection. Namely we have to keep all charge 2 Abelian 
Wilson loops in the effective action written by Abelian 
link fields to reproduce the screening of charge 2. Need- 
less to say, such an Abelian effective action is useless. 
The same problem is more serious in the real full QCD, 



since a fundamental charge also is screened. The authors 
in Ref. U suggest that the relevant quantity in the con- 
finement mechanism is not the Abelian monopoles but 
the Z(N) center- vortices which can explain the screen- 
ing problem |J. 

The aim of this note is to show how the screening and 
confinement problem is solved qualitatively in the frame- 
work of [/(l)^" 1 Abelian dynamics. We find the effect 
of dynamical charged particles can be described in terms 
of integer electric currents on the lattice. 

Let us first review the case without a dynamical char- 
ged particle. For simplicity, we consider only SU (2) QCD 
— ► U(l) and use the differential forms on the lattice |l0[ . 

Consider a modified compact QED action of the form 



z x = / ve 



-{d8+2Trn,AD(d8+2irn))+i(QJ,8) 
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where the operator D is a general differential operator, Q 
is a charge of an external source and the current J takes 
±1 along the Wilson loop. A is the Laplacian on the 
lattice. Following Ref. [|| one can reduce the partition 
function (m) to that of the following monopole action 



Zn 



E 



-(* k,D* k)+2-!ri(* k* N)—r 2 Q 2 (J,(A 2 D)- 1 J) 



(2) 



«fce^(* ci ) 
S*k=0 



where *N = A^^-S* S and S is an open surface spanned on 
the Wilson loop, i.e., 8S = QJ. Here and below the con- 
stant pre-factors are ignored. Numerically it is found Q] 
that the Abelian effective monopole action derived from 
SU(2) QCD can be described well in terms of two-point 
interactions alone in the long-range region as in (Q) . 

It is also possible to transform the monopole partition 
function (||) into the string representation as follows [||: 



Z 3 = 
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where the summation goes over the surfaces spanned (\Q\ 
times) on the external current J. 
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In the infrared limit the operator D is approximated 
well by Coulomb+self+ nearest neighbor terms Q|, that 
is, D = /?A _1 + a + 7A, where (3, a and 7 are renormal- 
ized coupling constants of the monopole action. They 
satisfy the relation (3 ^> a, 7. Then the Wilson loop can 
be estimated from Eq.(||): 



(W Q (R,T)) = const, e 



-« |Q| RT+ — 



(4) 



where i?T is the area of the minimal surface spanned on 
the contour J. There are \Q\ such surfaces which must 
be parallel to each other to maximize the contribu- 
tion for \Q\ 5s 2. This explains linearity in \Q\ in Eq. (^). 
When /? > 0,7, the string tension is evaluated approxi- 
mately as k = TT 2 /f3. 

Now let us introduce a dynamical charged particle. A 
charged vector field corresponding to gluon fields and 
a fermion field corresponding to quark fields are rather 
complicated to deal with. Hence we consider first a 
charged scalar Higgs field in the London limit as a sim- 
ple example. The radial part p x of the scalar Higgs field 
is frozen and the dynamical variable is the 
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compact Higgs phase $ € [—if, 71") which carries the elec- 
tric charge q. The Higgs field action is written in the 
Villain representation: 
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where G is a local operator. 

The integration of the phase of the Higgs field # can 
be represented as a weighted sum of the Wilson loop over 
the closed charged currents following Ref . [Q : 



/TT /*00 
m vf 
-7T J — OO J f- rjp f „ 

= fm £ 

J — 7T ■ ^ m / 



-(F^AGy 1 F)+i(F,di}+q9+2-nl) 



E < 

jez;(ei) 
5j=0 



= -(i,(4G)- 1 j)+ l (M,e) 



(6) 
(7) 



where we consecutively introduce the Gaussian integra- 
tion over the non-compact link variable F, apply the 
Poisson summation formula, e 2 ' Ki ( l > F ' = £\ S(j — F), 
integrate over the fields F and The last integration 
gives the closeness constraint Sj = 0. 

Hence the Villain type compact QED with the charged 
scalar field is written as 
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This expression can be reduced further to the monopole- 
electric and the string-electric current models: 
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where 



A X (5*s and s is a surface spanned q times 



on the dynamical current j, i.e., Ss = qj. 

FromEq.(|), we can see how the screening-confinement 
problem is solved. If Q/q £ the static potential can 
not be screened completely. For example, in the case of 
Q = I and q = 2, the leading string tension is equivalent 
to the string tension without presence of the dynamical 
charges as in (|J). In this case, j = in the sum over 
j gives the leading term while other terms coming from 
non-zero j show stronger damping. 

On the other hand, if Q/q = N e "Z (as in the case of 
Q = 1 and q = 1) then the expectation value of the RxT 
Wilson loop is expanded as a perimeter term (oc (R + T) ) 
given by the action of the dynamical screening current 
with QJ + qj = and the area-law term: 



(10) 
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These terms correspond to the screening currents, respec- 
tively, j = — N J and j = (— JV+1) J. The string breaking 
is seen at large distances. It is important to note that the 
area-law behavior is observed in the intermediate region 
even when the screening occurs. 

One can easily show why the screening is better ob- 
served in the case of the Polyakov loop correlator rather 
than in the case of the Wilson loop quantum average as 
found recently 0. Consider two Polyakov loops sepa- 
rated by a distance R corresponding to a pair of static 
quark and anti-quark in the case of periodic boundary 
condition in the time direction. The two leading terms 
in the average are: 

(P Q (0)P Q (R)) ^ d e- mN2T +d 1 e-^ RT - m ^ 2T +- ■ ■ . 

A comparison of this expression with Eq. ( |l0| ) gives that 
the area-law terms are the same for both averages. How- 
ever, the perimeter terms are different: for the Wilson 
loop the perimeter term contains additional suppression 
factor ~ e -™(2JV-i)fl w ith respect to the area term. Thus 
at sufficiently large separations between the sources the 
perimeter term in the Wilson loop average may not be 
found numerically. As a result, the string breaking may 
not be observed in the Wilson loop even if the breaking 
can be seen in the Polyakov loop correlator. 
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To derive numerically such an effective action (||) 
in terms of monopole and electric currents is very in- 
teresting. Integer monopoles can be defined following 
DeGrand-Toussaint[|l3| while definition of the integer 
electric currents in terms of original Abelian link fields 
seems to be impossible. Nevertheless, a given effective 
matter action as in (^), the action in terms of integer 
electric currents (^|) can also be derived. 

Now we consider the case of dynamical charged link 
variables which correspond to off-diagonal gluon fields 
after the Abelian projection of QCD. The SU(2) link 
variable can be parameterized as U Xlii = c x ^u x ^ where 
the matrix c Xi{i corresponds to the off-diagonal gluon 
represents the diagonal photon contribution: 



while u Xlli 

r ll _ „22 
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e l8 ^» and u x 2 ^ = u x \ = 0. The link fields 



x.^i x,fi 

are restricted, -q G [0, n/2], 8 7 ip e 



-7T,7rj 



The meaning of these fields can be understood after 
applying the U(l) gauge transformations, which are al- 
lowed in an Abelian projection. The field is the Abelian 
gauge field and p> is the vector field (phase of the off- 
diagonal gluon field) with the electric charge two, re- 
spectively: 6 X ^ -> 9 X ^ + £ x ~ Zx+fi, <Px,n -* fx.,^ + 2£ x . 

In the continuum limit these transformations reduce to 
0n(x) -> 0^{x) - d^(x) and ip^x) -» ip^x) + 2£(x), 
respectively. The variable r\ x ^ is U[X) gauge invariant. 

The SU(2) plaquette action Trllp in terms of the angles 
?/, 9 and p can be written as a sum of three parts [55 . 
The first part is proportional to the compact QED ac- 
tion, S e = -(3su(2) cosd6» Yli=i C0S Vi, where (3 S u(2) is 
the SU(2) gauge coupling and the subscripts 1, ...,4 cor- 
respond to the standard notations of the plaquette links: 
1 — » {x,x + //},..., 4 — > {x,x + £*}. The second part, 
S 8ip , contains 6 terms representing the interaction of the 
matter field p> with the gauge field 6. All these terms 
are proportional to the products cos 771 cos 772 sin 773 sin 774 
and the like. The last part of the action is given by 
= -f3 SU (2) cosd((p + 6) fj i=1 sinr^, where dip = 
Pi - P>2 + Pz - Pi- 

The simplest way to demonstrate the adjoint string 
breaking is to use a mean-field approximation for the 
field 77 (for our purposes the fluctuations of this field 
are not essential). We set cos 77^ — * {cosri x ^ L ) = c 
and sin 772; M — > (sin r/ Xtf j_) — s, respectively. Then the 
self-interactions of the gauge and matter fields can be 
written in a short form, S = — fisu{2)C A cosd6? and 
= — /?5(7(2)S 4 cosd(( / 9 + 8), respectively. 

The role of the actions S 6v and S v can be elucidated by 
noticing that they basically contains the following com- 
binations of the fields: self-interaction of the gauge field, 
d8, the interaction of the gauge field with the off-diagonal 



field, H. 



x,i±v 



2dx,ii+Px,i^~Px+fi,u and self-interaction of 
the off-diagonal field, C x ^ v = p x ^ — p Xl v Thus from the 
bare Wilson gauge action we get the following mean-field 



types of interactions involving the off-diagonal fields: 



qOtp _|_ Qif 
^X.^LU ~ ^X^LV 



cos_ff r 



cos H r 



cos C x 



The tensor is not antisymmetric contrary to C M „. 
Therefore we need two terms with H at each plaquette. 

A general effective action of charged link variables can 
be written in terms of two periodic actions: 



S 9 « = F H {H) + F H (H') + F C {C) , 



(11) 



where the form H = H Xjtlu {H 1 = H x v ^) corresponds to 
the positively (negatively) oriented plaquette interactions 
of the matter and gauge fields. Now let us consider the 
adjoint string breaking in the SU(2) gluodynamics using 
the following partition function: 



Z 7 = V6 Vipe 



(12) 



The screening of the external charge-2 (adjoint) Wilson 
loops is possible essentially due to the presence of the 
off-diagonal field action S 6v . To show this we apply the 
duality transformation with respect to this action: 



V J H T H T C P ^ H ^ 



)+i(H' ,n 2 )+i(C,n 3 ) 



i=l,2,3 



where the weights /„ . correspond to the terms in Eq. <h 



The integration of Eq. (12) over p> gives us the sum 
over the particle trajectories with the charge 2: 



Vp e 



E 

Sj=0 



w(j) 



(13) 



where w(j) is the weight corresponding to the ensemble 
of the Wilson loops {j}. The integration over the field 
ip provide us with constraints on the dynamical integer 
trajectories j: the trajectories must be closed, Sj = 0, 
due to gauge invariance of both parts of Eq. (|l3|). 

The final step is to substitute the loop expansion {Jlc 
into the partition function ( |l2| ) and to apply the argu- 
ments identical to those in the charge Q — 2 scalar par- 
ticle considered earlier. Thus we get the screening of the 
charge-2 external charges (gluons) and, as a result, the 
breaking of the adjoint string. 

Let us next discuss the case of charged fermions corre- 
sponding to quark fields. The fermion determinant can 
be rewritten as a sum of self-avoiding Wilson loops com- 
posed of charge q = 1 particle [l6|. Hence we get 
an action similar to (|^) and the screening-confinement 
problem could be solved similarly. 

But generally an effective fermionic action becomes 
more complicated and moreover it is desirable to adopt a 
form of action which makes numerical calculations pos- 
sible. We employ a Polynomial Hybrid Monte Carlo 
Method lO. 
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Consider a polynomial -P/v , [z] with an even degree where wf is a weight for the electric particle trajectories. 
Npoiy which approaches z~ l as N po i y increases. We get 



det[P N ,[D}] = det 



N„, 



[D] (14) 



where Tn t [z] is a polynomial with a degree N po i y /2. 
Hence the fermionic determinant for Nf odd-number of 
flavors is written using a pseudo-fermion representation: 

,v _ r^rn D .. rnlll^ (15) 



(det[D])' 



[det[DP Npoly [D]}]^ 



exp 



Tk^[D]T Npol \D]<t> 



Since T^ volv [z] is a polynomial of degree N po i y /2, then 



Vpoly L 



[D] =a + a x M + a*M t 



(16) 



where M is a hopping parameter term in D. Since the 
higher-order terms are irrelevant in the screening prob- 
lem, we consider only the linear term with respect to M. 
Here we adopt the Wilson fermion action for M coupled 
to a U(l) compact gauge field as an example: 



M = K 



4 

E 



{(1-7/.) 







+ (1 + 7m) 











(17) 



We separate the U(l) variant phase factor from the U(l) 
invariant term (f> a (x) of the pseudo-fermion field <j> a (x) = 
4>a{x) exp(ix XiQ T3) (a =Dirac indices). Regarding the 
field <j> as a constant we get the following mean-field types 
of interactions from the Wilson-fermion action: 



^ cos(A m Xx,q + 9 x ,m) > 



a,fi,x 

E C0S (Xx+ M ,/3 - Xx,a + O x ^)%0(ji) , 
a,/3,^i,x 

where 7(/i) are realized Euclidean gamma matrices, i.e., 
7(/i) = i7 M for /i = 1, 3 and j(fi) = 7 M for fi = 2, 4. 

As an effective fermionic action, it is natural to assume 
the following periodic pseudo-fermion action: 

p -S f {e,x) _ „-F(A tlXx , Q +e x , fl )+G(xx +M ,^-Xx, Q +e :c , fI )7a3(M) 



x e 



E E 

i Q . M ( a: )G2ic«3,K( a: )e2 



FUa^GUa^) (18) 



x e 



y'a/3,f.( :c )(Xa!+M,3-X ! c,a+0a ! , / i)7a/3(/i) 



where F and G are periodic functions. In Eq. ( |18| ) we 
have used the Fourier (duality) transformation. The in- 



teger current j M (x) = S«ic 



is coupled to the field 6 and carries unit electric charge. 
Due to the gauge invariance of the fermion determinant 
the integration of Eq. ( |l8| ) over the phases x gives us a 
conservation constraint for the currents, 5j = 0. Thus 



V X e 



-s f (e,x) 



(19) 



With the gauge field action (|I|), we can see that the 
screening and the string breaking in full QCD also is suc- 
cessfully explained as done in the case of charged scalar 
particle in the framework of U(l) theory. It is again 
stressed that the area-law behavior is seen in the inter- 
mediate regions even in full QCD. Such an area-law be- 
havior can be explained by the dual Meissner effect due 
to the monopolc condensation. Confinement is important 
also in full QCD. 
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